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Abstract We systematically compare filtering methods used to extract topological excitations (like in- 
stantons, calorons, monopoles and vortices) from lattice gauge configurations, namely APE-smearing and 
spectral decompositions based on lattice Dirac and Laplace operators. Each of these techniques introduces 
ambiguities, which can invalidate the interpretation of the results. We show, however, that all these meth- 
ods, when handled with care, reveal very similar topological structures. Hence, these common structures 
are free of ambiguities and faithfully represent infrared degrees of freedom in the QCD vacuum. As an 
application we discuss an interesting power-law for the clusters of filtered topological charge. 
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1 Introduction 

Ever since the advent of quantum chromodynamics (QCD) 
its infrared properties have been of primary interest. One 
of the most important phenomena in this regime is the 
confinement of quarks. As a typical nonperturbative ef- 
fect it still calls for a derivation from first principles, even 
for pure gauge theory (where it reveals itself as, e.g., an 
area law for the Wilson loop). 

Remarkably, most of the most popular nonperturba- 
tive approaches to the QCD vacuum have been around 
for quite some time now. Topological excitations like in- 
stantons, calorons, magnetic monopoles and vortices have 
been used for semiclassical and condensed matter-inspired 
models since the 70's pQ. In spite of the successes of these 
models, the question of their relevance for confinement is 
the subject of debate. 

Lattice gauge theory (LGT) is of roughly the same age 
and till today the only nonperturbative and systematically 
improvable regulator of QCD. Thus it has the potential to 
lend support to and to decide between the various models. 

However, the QCD vacuum on the lattice as seen by 
naive gluonic observables has been found to be dominated 
by ultraviolet (size 0(a)) fluctuations and therefore it is 
difficult to make contact to continuum models. To deal 
with this problem, various smoothing procedures, filtering 
out the UV 'noise', have been developed and applied to 
lattice configurations. It has been objected that physical 
properties could be lost and unphysical artefacts may be 
generated by these filtering methods. Both effects would 
strongly spoil the conclusions, for instance w.r.t. the ex- 



tracted density of the building blocks, drawn from any 
such study. 

In this work we systematically investigate the most 
common and a priori quite different procedures to filter 
lattice configurations, namely smearing and the modern 
methods based on the eigenmodes of lattice Dirac and 
Laplace operators. 

We find a surprisingly strong agreement of the topo- 
logical content of configurations in thermal equilibrium 
seen through the different methods. We show how the pa- 
rameters of the latter can be adapted systematically. This 
knowledge forms a very important prerequisite to uniquely 
identify the structure of the QCD vacuum. 



2 Filtering methods 

Smearing and the related method of cooling have often 
been used to improve the signal in observables or to ob- 
tain smoother link variables as input for lattice operators. 
For definiteness here we will use a 4D version of APE- 
smearing, that has been argued to be equivalent to RG 
cycling [2] . It is an iterative procedure, where links are re- 
placed by a weighted average of the links and the staples 
U"(x) = U u {x)U^{x + £>)Ul(x + p) surrounding it: 

U^x)^V[aU,(x) +1 Y, U ^ x )]- (!) 

Here V denotes the projection onto the gauge group (for 
577(2) just a rescaling of the matrix by a scalar). We 
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Figure 1. Effect of different filtering methods on the topological density for a particular Q — 1 configuration in a fixed lattice 
plane. On top we show the original topological density. In the second and third row one sees the effects of smearing (first 
column) and Laplace filtering (second column) as well as the topological density in terms of Dirac eigenmodes (last column). 
At the bottom the profile of the chiral zero mode is shown for the combined filtering methods. (The numbers in brackets give 
the heights of the corresponding maxima. 



choose a = 0.55 and 7 = 0.075, following [2]. Cooling 
is obtained by ignoring the old link (a = 0) and is known 
to drive the configurations towards classical solutions. We 
update one link at a time, but as the weight of the cen- 
tral link is quite high, our technique qualifies as smearing 
rather than cooling. 



The Dirac filtering method relies on the representation 
of gluonic observables through eigenmode expansions of 
lattice Dirac operators (see also (4|). In the following we 
will investigate the topological charge density [5] in terms 
of eigenmodes of a Ginsparg- Wilson type Dirac operator 
D: 



A more recent idea for filtering is to use eigenmodes 
of lattice Dirac operators. Generally speaking, their func- 
tion as filters is based on the argument that low-lying 
eigenmodes tend to be smooth, see Fig. Q] for an example. 
The positions revealed by the lowest-lying modes are ex- 
pected to be correlated with the location of the relevant 
gluonic IR excitations, in particular of topological objects 
[3]. Whether peaks show up at the same locations when 
gluonic filtering methods are applied, is an important con- 
sistency check for the latter. 



1 A 

- tr l5 (-D x , x -l) = J2^- l)^(*)7^n(a0 (2) 



which is exact for N — Vol • 4jV c and can be truncated 
for filtering purposes at low N [6]. The total topological 
charge Q = q(x) is an integer given by the contribu- 
tions of chiral zero modes, whereas the non-zero modes 
modify the local distribution q(x) only. 
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The third filtering method we use is the Laplace filter 
that represents link variables through a spectral sum [7] 

JV 

U^x) = V^{l&-2\ n )<t> n {x)®<t>t{x + (i). (3) 

71=1 

Here, <j) n are the eigenmodes of the gauge covariant Laplace 
operator —D 2 [U] with eigenvalues A„, of which again only 
the lowest N are taken into account. In the limit N = 
Vol • N c the original links would be reproduced. Numeri- 
cally this filter is cheaper than Dirac operators. 

The strength of the filters can be controlled. More iter- 
ations of smearing or fewer eigenmodes result in a stronger 
filtering. We also remark that smearing is a strictly local 
procedure, while working with eigenmodes introduces an 
intrinsic nonlocality, the consequences of which are not 
fully understood yet. 
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10 or more modes we find again one chiral zero mode con- 
firming Q = 1. 

The example discussed so far is typical for the whole 
ensemble (including all Q-values) and helps one to visu- 
alize the impressive similarity of smearing, Laplace filter- 
ing and Dirac eigenmodes, which in the following will be 
quantified. 



3 Results 

We have generated 295 independent quenched SU(2) con- 
figurations on a 16 4 lattice with tree-level Liischer-Weisz 
action at (3 = 1.95 (the lattice spacing is a = 0.075(1) fm). 
We use chirally improved fermions [5] , an approximate so- 
lution of the Ginsparg- Wilson relation, which reveal chiral 
properties well enough without the big computational de- 
mands of, e.g., overlap fermions. 

The main physical observable for our comparative study 
is the topological charge density. For smearing and the 
Laplace method, which both provide filtered links, we use 
the gluonic definition q(x) = £\ v trF llv F lxv /l&'K 2 with 
an improved field strength F^ u constructed from lxl, 
2x2 and 3x3 Wilson loops [9] . The Dirac eigenmodes 
yield q(x) directly via Eq. @. 

In Fig. Q] we show the result of the three methods real- 
izing two levels of filtering (mild and strong, see later) on 
a thermalized Q = 1 configuration in a fixed lattice plane. 
It shows an excellent agreement of the 'hot spots', i.e., 
lattice locations with large local topological charge, visi- 
ble through the various methods. For mild filtering also 
less pronounced structures appear, which agree between 
the methods (and are washed out by strong filtering). 

Moreover, the most pronounced structure corresponds 
to a maximum in the profile of the chiral zero mode. From 
this consistency we can conclude qualitatively that differ- 
ences between the methods are small. 

The configuration used in Fig. Q] has unit topological 
charge, following from the existence of a single chiral zero 
mode. The following table shows in as far smearing and 
Laplace filtering recover this fact: 

The main conclusion here is that the three methods point 
to the same total charge Q = 1. Furthermore, Q evaluated 
after smearing approaches unity rather quickly and then 
becomes stable. It also keeps the chiral zero mode (marked 
with asterisks * in the table). The Laplace-filtered config- 
urations, on the other hand, are not arbitrarily smooth, 
such that the gluonic topological charge measured on them 
is 1 within some margin; on the configurations filtered with 




Figure 2. Auto-correlators XAa(t) (broken lines) and cross- 
correlators xab{t) (full line) of smearing (S), Laplace (L) and 
Dirac (D) filtering for the configuration of Fig. Q] at mild filter- 
ing. The unit on the vertical axis is 10 -7 . 



In particular we try to set the parameters of the fil- 
tering methods such that the results match as well as 
possible. For that purpose we consider correlators of the 
topological charge density (with its average q = Q/Vol 
subtracted): 



Xab{t) 



T, x J<1a(x) - q A )(q B {y) - q B )S{\x - y\ - r) 



J2 x , v s (\ x -y\ - r ) 



(4) 

depending on the four-dimensional distance r. A and B 
stand for the filtering methods under consideration (in- 
cluding their parameters). 

As Fig. [2] shows, the auto-correlators Xaa{t) have a 
positive profile over a few lattice spacings followed by a 
slightly negative tail, and the cross-correlators xab (0 fall 
in between them. The ratio of the latter to the geomet- 
ric mean, Sab = Xab( q )/(xaa(0)xbb(0)), being close to 
1 signals local agreement between the topological charge 
landscapes of methods A and B. 
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Figure 3. Level curves Sab = 0.1, 0.2, . . . (inwards) for A and 
B being smearing sweeps (S), Laplacian (L) and Dirac (D) 
eigenmodes, averaged over 10 configurations. The parameter 
sets we use at mild and strong filtering are indicated by a full 
triangle and a circle, respectively. 



Fig. shows Sab for pairwise comparisons of smear- 
ing, Laplacian and Dirac filtering]. In all pairs of methods 
Sab reveals a 'ridge' on which the parameters match best. 
From mild filtering (lower left corneJJ) to stronger filter- 
ing (upper right corner) the methods deviate from each 
other and the height of the ridge decreases. Two optimal 
parameter sets chosen from that figure (and indicated in 
it) are used in due course: 



Investigations that need 50 Dirac modes are done on an 
ensemble of 10 configurations for computational reasons. 

2 From the definition it follows that in the limit of no smear- 
ing sweeps vs. Vol -Nc Laplacian modes one has Sab = 1. 



mild filtering 
strong filtering 



smearing 
10 sweeps 
20 sweeps 



Laplace 
80 modes 
20 modes 



Dirac 
50 modes 
8 modes 



These matchings have also been confirmed by comparing 
smeared and Laplace-filtered links minimizing 



UZ(x))Hu£(x)-U?(x)) 



For these two optimized parameter sets we summarize in 
the following table some additional measurements, which 
characterize the deviation of the smeared and Laplace- 
filtcrcd ensemble from the orginal Monte Carlo one (295 
configurations) . The observables are the percentage of con- 
figurations for which the topological charge Q coincides 
with the number of zero modes (called Qd) and the de- 
crease of the action S: 



10 sweeps 
80 modes 



Qd \ < 0.5 

89% 

76% 



S I S'orig 

0.026 
0.036 



20 sweeps 
20 modes 



85% 
69% 



0.009 
0.017 



In reducing the noise, the action has been reduced by ap- 
proximately two orders of magnitude, while the string ten- 
sion is unchanged within errors. 

With the parameter sets at hand we further quantify 
the local agreement of the methods by looking at clusters 
of filtered topological charge. For each method A, cuts in 
the topological charge density are adjusted such that the 
sets Xa of points above the cut have the same volume 
fraction / = volpC^/Vol. Then we compare the volume 
in the overlap (with same sign of topological charge) vs. 
the union of pairs Xa and Xb to obtain the relative point 
overlap (RPO) sab' 



sab 



E 



x<=x A nx B 
qA(x)qB(x)>0 



w E i- 

x&X a UXb 



(5) 



In Fig. |4] upper panel the RPOs are plotted against the 
volume fraction /. The main result to emphasize here is 
that the pairwise point overlaps are large, typically 50% 
to 60%, and constant over a wide range of /, which implies 
that also the shapes of the topological lumps agree. 

Based on our analysis we are able to discard ambi- 
guities of the individual methods. When we analyze only 
points that are common to all three methods, the number 
of clusters is reduced accordingly, see Fig. |4] lower panel. 

Intermediate values of / arc suited for a cluster anal- 
ysis as they contain enough statistics and avoid cluster 
mergings. Interestingly, in that regime we find for the com- 
mon structures a pronounced power-law for the number of 
clusters as a function of f. In order to better interpret this 
behavior, we will discuss in the next section gross features 
of a class of models including the dilute instanton gas that 
predict such a power-law. 
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Figure 4. RPOs averaged over 10 configurations (left) and the number of clusters for the configuration of Fig. \T\ (right, 
logarithmically, including a power-law fit), both as a function of the volume fraction / at mild filtering. 



4 A model for the power-law 

We are now considering models of topological lumps in 
the continuum. The latter are assumed to be dilute and 
described by an arbitrary shape function F with size pa- 
rameter p: 

q(r) = F(r/p) p- s . (6) 

When applying cuts to the topological charge density, it 
is useful to characterize the lumps by their maximum 

q Q = q(0) = F(0)p- s . (7) 

Let us assume a power-law distribution of size parame- 
ters d(p) ~ pP , which can be translated into a power-law 
distribution of maxima 

d(<fo)~Sb" , <* = (/3 + 5 + 1) /S. (8) 

Then the number of clusters visible above a cut q is simply 

N ciust (q) = / dq d(q ) ~ q 1 ^ . (9) 

Jq 

The next step is to compute the total volume of these 
clusters. We assume that r is a <i-dimensional radius. Then 
for a single cluster (with maximum qo) we obtain the fol- 
lowing volume above q: 

V(q ,q) = j d^Q J" ' dr r*- 1 e(q(r) - q) (10) 
~ p d F-\qp s ) ~ q^ d/5 F-\F{Q)q/q Q ) . 

Hence the total volume or equivalently the total number 
of lattice points with topological charge density above the 
cut q can be calculated to 

iV P oints(<z) = ^ dq d(q )V(q, q ) ~ q 1 ^' 5 . (11) 
Jq 



Together with I|8I9[) this predicts the following exponent 

dlogjV c i ust (g) = 1 
? " dlog^ points (g) l + d/(/3+l) ' 1 1 

independently of the actual profile F. 

For the naive dilute instanton ensemble the various 
parameters are d = 8 = 4 and (3 = lliV c /3 — 5 = 7/3, 
yielding £ = 5/11 (the precise way of cutting off the in- 
stanton size distribution at large p will not matter since we 
are concentrating on hot spots of large topological charge 
density, i.e. small size). 

Tests with toy models of lumps show that the detection 
of the power-law and its coefficient £ hardly depends on 
finite volume or finite lattice spacing, as long as the lumps 
are dilute. However, a high density of topological objects 
modifies £. In the extreme case of the topological charge 
density being pure noise one obtains £ = 1, as every cluster 
comes with just one point. 

The measured exponent for the topological charge den- 
sity indeed is close to 1 for very mild or no filtering. As 
Fig. [5] shows, £ is lowered when removing noise by filter- 
ing and tends to a plateau for moderate filter parameters. 
At strong filtering the cluster statistics becomes unreliable 
since the topological landscape is dominated by very few 
clusters. 

In order to avoid ambiguities of the different methods, 
we have based our analysis on clusters common to smear- 
ing and Laplacian modes (neglecting the numerically ex- 
pensive Dirac eigenmodes). In the upper panel of Fig. [5] 
we depict £ along the ridge of optimally matching param- 
eters, whereas in the lower panel of that figure we include 
other (non-optimal) pairs of parameters. 

While the statistical errors of the fits are small, the sys- 
tematic errors from choosing the /-interval for the power- 
law are hard to quantify. Estimating the latter we obtain 
the following value of the power-law exponent 

f = 0.59(5) (13) 

for the filtered lumps of topological charge density. 
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Figure 5. The exponent £ of common topological clusters as a function of smearing sweeps and Laplacian modes. The best 
matching parameters (cf. Fig. [3] right) are used in the left plot and marked by a grey band in the right plot. 



This result certainly excludes the simplest dilute gas 
of instantons (for which £ = 0.45 see above, included in 
Fig. [5] top) as primary topological structures. This is not 
surprising, as instanton ensembles are expected to be non- 
dilute. For other scenarios the expected behavior depends 
crucially on the specific model assumptions. Therefore, we 
leave it to the proponents of any such model to check 
whether it is in agreement with our findings. 

5 Conclusions 

We have presented a systematic comparison of filtering 
methods on the lattice, namely the gluonic smearing up- 
date and truncated expansions based on Laplacian and 
Dirac eigenmodes. Our finding, that these filtering meth- 
ods agree surprisingly well, allows us to identify structures 
of interest nearly free of ambiguities. 

As a general rule, smearing - which is the cheapest 
method - represents these objects reasonably well. In or- 
der to control its effects, however, we consider it essential 
to compare several filtering methods. To this end we have 
shown how to match the parameters of the methods opti- 
mally. 

As a first application we have measured a power-law 
for the number of filtered clusters and presented a class of 
models in which the corresponding exponent is a function 
of the size distribution coefficient and the dimensionality. 
We have found that the measured exponent cannot be 
interpreted in terms of a dilute instanton gas. 

This work has been supported by DFG (Forschergruppe 
'Gitter-Hadronen-Phanomenologie') and BMBF. 
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